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Abstract 


The  generalized  lattice  point  (GLP)  problem  provides  a  formulation  that 
accommodates  a  variety  of  discrete  alternative  problems.  In  this  paper  we  show 
how  to  substantially  strengthen  the  convexity  cuts  for  the  GLP  problem.  The 
new  cuts  are  based  on  the  identification  of  "synthesized"  lattice  point  condi¬ 
tions  to  replace  those  that  ordinarily  define  the  cut.  The  synthesized  condi¬ 
tions  give  an  alternative  set  of  hyperplanes  that  enlarge  the  convex  set,  thus 
allowing  the  cut  to  be  shifted  deeper  into  the  solution  space.  A  convenient 
feature  of  the  strengthened  cuts  is  the  existence  of  linking  relationships  by 
which  they  may  be  constructively  generated  from  the  original  cuts.  Geometric 
examples  are  given  in  the  last  section  to  show  how  the  new  cuts  improve  upon  those 
previously  proposed  for  the  GLP  problem. 


1.  INTRODUCTION 


The  generalized  lattice  po4nt  (GLP)  problem  provides  a  formulation  that 
accommodates  a  variety  of  discrete  alternative  problems.  Significantly,  it 
accomplishes  this  without  the  necessity  of  introducing  integer  variables  and 
their  associated  constraints.  The  advantage  of  such  a  formulation  lies  not 
only  in  the  reduced  problem  size  it  affords,  but  also  in  the  fact  that  solu¬ 
tion  methods  for  the  GLP  problem  do  not  have  to  wrestle  with  "extraneous  struc¬ 
ture"  created  by  standard  tricks  for  formulating  such  problems  as  integer  programs. 

A  class  of  cutting  planes  that  apply  directly  to  the  GLP  problem  formulation  utili¬ 
zing  the  convexity  cut  framework  was  first  developed  in  [2].  Extensions  of  these 
to  the  related  structures  of  "facet  problems,"  have  also  been  given  in  [3l> 

In  this  paper  we  show  how  to  substantially  strengthen  the  convexity  cuts 
for  the  GLP  problem.  The  new  cuts  are  b*aed  on  the  identification  of  "synthe¬ 
sized"  lattice  point  conditions  to  replace  those  that  ordinarily  define  the  cut. 

The  synthesized  conditions  give  an  alternative  set  of  hyperplanes  that  enlarge 
the  circumscribing  convex  set,  thus  allowing  the  cut  to  be  shifted  deeper  into 
the  solution  space. 

A  convenient  feature  of  the  strengthened  cuts  is  the  existence  of  linking 
relationships  by  which  they  may  be  constructively  generated  from  the  original  cuts. 
Geometric  examples  are  given  in  the  last  section  to  show  how  the  new  cuts  improve 
upon  those  previously  proposed  for  the  GLP  problem. 

2.  THE  GLP  PROBLEM 

The  GLP  problem  is  an  ordinary  linear  program  subject  to  an  additional 
"lattice  point  condition."  This  condition  stipulates  that  all  admissible  eola¬ 
tions  must  lie  on  an  (m-q)  dimensional  face  of  some  convex  polytope  defined  by 
a  subset  of  the  constraints  of  the  linear  program,  where  m  is  the  full  dimension 
of  the  polytope  and  q  is  a  constant  between  0  and  m.  Stated  in  the  language  of  basic 
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solutions,  feasibility  for  the  GLP  problem  requires  that  at  least  q  slack  variables 
associated  with  a  specific  set  of  constraints  be  excluded  from  the  linear  programming 
basis.  Thus,  the  problem  may  be  expressed 
Maximise  cx 
Subject  to  Ax  <  b 
Dx  <  d 

and  at  least  q  slack  variables  for  Dx  <  d  cure  nonbasic. 

The  vector  of  slack  variables  for  Dx  <  d  may  be  written  u  ■  d  -  Dx,  and  hence 
the  lattice  point  condition  requires  q  of  the  components  u^  of  u  to  be  nonbasic. 
(Nannegativity  conditions  such  as  x  >  0  may  be  assumed  to  be  absorbed  into  the 
inequalities  Ax  <  b  and  Dx  <  d.) 

A  variety  of  problems  can  be  expressed  in  this  form,  including  the  0-1  mixed 
integer  programming  problem  (see  [2,3])  and  all  problems  that  can  be  formulated  in 
terms  of  the  latter.  For  some  of  these  problems,  however,  the  GLP  formulation  is 
more  direct,  suggesting  that  an  approach  that  is  specially  designed  for  this  pro¬ 
blem  may  have  advantages  over  methods  designed  for  other  formulations.  In  particular, 
it  is  shown  in  [3]  that  0-1  mixed  integer  formulations  of  certain  problems  give 
rise  to  a  "structural  redundancy"  --  requiring  an  undetermined  number  of  cuts  to 
move  —  which  is  avoided  by  the  more  direct  formulation. 

A  consequence  of  the  requirement  that  q  of  variables  be  nonbasic  is  that  q 
of  these  variables  receive  a  0  value.  For  many  problems  (including  the  0-1  pro¬ 
blem)  these  conditions  are  equivalent,  and  we  will  focus  here  on  this  simple  version 
of  the  lattice  point  condition.  Our  results  apply  as  well  to  the  "nonbasic"  version 
by  the  procedures  indicated  in  [2]. 

3.  3LP  CUTS 

A  central  result  for  the  Gif  problem  states  that  one  can  select  any  set  of 


m  -  q  +  1  of  the  variables  and  determine  the  convexity  cut  relative  to  the  convex 
region  given  by  >  0  for  these  selected  variables.  Here  m  is  the  total  number  of 
u.  variables.  (The  rule  for  generating  the  cut  is  to  extend  each  edge  associated 
with  the  current  LP  tableau  until  it  intersects  the  boundary  of  the  convex  region, 
and  then  pass  the  cut  hyperplane  through  the  endpoints  of  these  edges.  See  [2]  for 
detailed  illustrations  of  how  this  may  be  done.  The  geometric  examples  of  Section  4 
also  provide  a  pictorial  analog.)  One  can  always  find  such  a  set  of  ^  variables 
all  of  which  are  positive  whenever  fewer  than  q  of  the  u^  are  zero.  The  positivity 
restriction  on  the  values  of  these  m  -  q  +  1  variables  assures  that  the  cut  is  not 
"degenerate"  relative  to  the  current  LP  extreme  point. 

To  determine  strengthened  cuts,  we  seek  new  sets  of  m  -  q  +  1  variables,  con¬ 
catenated  out  of  the  original  components  of  u,  which  have  the  requisite  properties 
to  define  suitable  convex  sets. 

The  result  that  identifies  these  new  variables  and  specifies  how  they  give  rise 
to  alternative  convexity  cuts  is  the  following. 

Theorem  1.  Partition  the  index  set  (l,...,m)  for  the  u^  variables  into  m  -  q  +  1 
nonempty  sets  Pj,  J  *  1, .  • .  ,m  -  q  +  1.  For  each  of  these  sets  Pj,  define  a  variable 

'j  -ifejVi 

where  the  parameters  are  chosen  so  that  each  Zj  is  positive  in  the  current 
basic  solution.  Then  a  legitimate  convexity  cut  for  the  GLF  problem  is  given  by  the 
convex  set  consisting  of  the  intersection  of  the  half  spaces  z^  >  0,  j»l,...,m  -  q  +  1. 

Proof:  Drawing  on  the  results  of  [1,2],  the  validity  of  the  theorem  is  easily  established. 
It  is  only  reqiired  to  demonstrate  that  at  least  one  of  the  variables  ij  must  be  0  in 
every  feasible  solution  to  the  GLP  problem.  Thus,  suppose  to  '-.he  contrary  that  there 
exists  a  feasible  solution  to  the  GLP  problem  in  which  all  of  the  s^  are  nonzero. 

This  implies  that  at  least  one  u^  in  each  set  P^  is  nonzero, 


and  hence  a  total  of  at  least  m  -  q  +  1  of  the  variables  are  nonzero,  fiat 
this  is  impossible  by  the  requirement  that  q  of  the  are  0  in  every  feasible 
solution. 

Note  that  this  foregoing  recult  specializes  to  the  earlier  result  for  the 
GLP  problem  by  imposing  the  restriction  that  only  one  X^  is  nonzero  for  each  P. 

(Thus  z j  >  0  corresponds  to  X^  >  0,  or  equivalently  u^  >  0;  for  each  of  m-q+1 
variables.)  We  new  indicate  how  the  theorem  can  be  used  to  generate  improved  cuts. 

Remark  1:  For  each  set  Pj  there  must  be  at  least  one  u^  (i  e  Pj)  that  is  nonzero 
in  the  current  basic  solution  in  order  for  the  Zj  all  to  be  positive.  Moreover, 
any  collection  of  m  -  q  +  1  positive  variables  can  be  taken  to  be  the  "starting" 

(or  representative)  variables  for  the  sets  P^.  Given  a  choice  for  these  starting 
variables,  the  remainder  of  the  sets  Pj  can  be  constructed  by  successively  assigning 
each  of  the  remaining  q  *  1  variables  to  them.  This  permits  the  cuts  from  Theorem  1 
to  be  generated  by  stepwise  modification  of  the  earlier  cuts  (i.e.,  their  correspond¬ 
ing  convex  regions). 

Remark  2:  Assume  the  basic  u^  variables  are  the  first  to  be  assigned  to  the  sets 
P  .  (The  starting  variable  for  each  set  will  of  course  always  be  basic.)  Then  the 
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constructive  process  of  Remark  1  can  be  simplified  for  each  of  the  nonbasic  uA  variables 
to  determine  the  "best"  set  Pj  to  which  it  should  be  assigned.  Moreover,  it  is 
possible  to  identify  the  deepest  cut  from  such  an  assignment  without  bothering  to 
specify  the  corresponding  X^  values.  In  particular,  suppose  that  "incomplete"  z^ 
variables  have  been  defined  relative  to  the  partially  constructed  F.  sets  (consisting 
of  the  basic  u^)  so  that  the  halfspaces  Zj  >  0  already  implicitly  define  a  cut.  To 
deepen  this  cut,  assign  each  nonbasic  u^  to  the  set  for  the  incomplete  z^  variable 
that  first  goes  to  0  as  a  result  of  extending  the  edge  corresponding  to  u^.  (This 
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Zj  variable  is  identified  by  the  same  ratio  test  used  to  identify  the  variable  to 
leave  the  basis  in  the  primal  simplex  method.)  The  best  value  will  produce  a 
new  Zj  variable  with  a  0  coefficient  in  its  "defining  equation"  (constructed  fran 
the  defining  equations  of  the  u^  variables  in  the  current  LP  tableau),  and  thus 
the  edge  corresponding  to  the  nonbasic  u^  variable  will  no  longer  be  blocked  by 
this  2^  variable.  Thus,  the  completed  cut  is  constructed  simply  by  permitting 
each  edge  corresponding  to  a  nonbasic  u^  variable  to  bypass  the  first  z^  that  blocks 
it.  (This  will  always  yield  a  deeper  cut  unless  more  than  one  z  qualifies  as  the 

u 

"first"  blocking  variable.) 


Remark  3:  Basic  ^  variables  may  be  added  to  the  Pj  sets  (given  a  selection  of 
m-q  +  1  starting  variables,  as  indicated  in  Remark  l)  by  similarly  considering 
incomplete  z^  variables  that  currently  block  particular  edges.  To  understand  hew 
this  may  be  accomplished,  note  that  a  linear  combination  of  an  incomplete  z^  with 
a  u^  creates  a  new  z^  variable  z^ '  whose  associated  hyperplane  z^ '  =  0  is  rotated 
through  the  intersection  of  the  hyperplanes  Zj  *  0  and  u^  ■  0.  Provided  the  hyper¬ 
plane  z.  ■  0  does  not  block  more  than  one  of  the  edges,  it  is  always  possible  to 
specify  a  Zj '  which  allows  the  edge  blocked  by  Zj  to  be  extended  an  additional  dis¬ 
tance  without  restricting  the  extensions  of  other  edges.  (This  is  a  simple  matter 
algebraically,  accomplished  by  examining  ratios.)  In  the  absence  of  "blocking  ties," 
such  a  procedure  results  in  a  strictly  deeper  cut.  (Linear  combinations  of  basic  u^ 
variables  may  alternatively  be  selected,  if  desired,  to  allow  certain  edges  to  be 
extended  still  more  deeply  at  the  expense  of  curtailing  other  edge  extensions.) 


Remark  4:  A  good  strategy  for  assigning  a  basic  u^  variable  to  an  incomplete  set  Pj 
is  to  select  a  variable  z^,  when  possible,  which  currently  blocks  a  given  edge  and 
which  would  continue  to  block  that  edge  if  u^  >  0  were  added  to  the  set  of  half 
spaces  defining  the  cut.  This  means  that  the  edge  intersects  u^  ■  0  after  lntersect- 
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lng  z  *  0,  tending  to  increase  the  depth  to  which  a  linear  combination  of  ^  and 
Zj  will  allow  the  given  edge  to  penetrate  before  interfering  with  the  extensions 
of  other  edges. 

More  complex  strategies  can  be  employed  for  generating  additionally  strengthened 
cuts  by  Theorem  1.  For  example,  basic  variables  whose  defining  equations  have  a 
number  of  oppositely  signed  coefficients  in  the  LP  tableau  may  be  used  to  "offset" 
each  other  in  the  context  of  Remark  3;  and  thus  may  be  paired  instead  of  considered 
independently  to  permit  deeper  edge  extensions.  Similarly,  the  effect  of  the  final 
strengthening  of  the  nonbasic  u^  variables  may  be  anticipated  when  determining  linear 
combinations  of  basic  variables.  Indeed,  having  determined  an  assignment  of  the 
variables  to  the  sets  Pj,  "optimal"  values  of  the  parameters  (e.g.,  those  which 
maximise  a  weighted  sum  of  the  depths  of  the  edge  extensions)  may  be  generated  by 
solving  a  simple  linear  program.  Such  refinements,  however,  are  not  required  to 
produce  cuts  that  are  stronger  than  the  original  GLP  cuts  since  the  procedures  of 
the  preceding  remarks  will  typically  suffice,  sometimes  dramatically  so. 

4.  ILLUSTRATIONS 

We  give  two  examples,  each  in  two  dimensions,  to  illustrate  the  new  cuts  and 
their  relationship  to  the  old  ones.  The  first  is 

Minimise  x^  +  xg 
Subject  to:  Xj^  +2xg  <  8 
x±  -2xg  <  2 

xr  x2  -  0 

where  the  lattice  point  conditions  require  at  least  2  of  u^,  Ug  and  u^  to  be  0, 
defining  ■  8  -  x^  -  2xg,  ug  ■  2  -  x1  +  2Xg  and  u^  ■  x^  (that  is  u^,  Ug  and  u^ 
are  the  slack  variables  for  the  first  three  Inequalities  of  the  problem). 


The  constraints  of  this  problem  are  depicted  in  Figure  1,  with  the  hyperplanes 
corresponding  to  u^  =  0,  ug  ■  0  and  ■  0  denoted  by  (l),  (2)  and  (3)  respectively. 

Thus  the  feasible  region  of  the  problem  excluding  the  lattice  point  conditions  is  the 
quadrilateral  bounded  by  (1),  (2),  (3)  and  the  x^  axis  (i.e. ,  the  hyperplane  xg  =  0). 

The  optimal  LP  vertex  is  the  origin.  The  GLP  cut  of  [2]  requires  selecting  2 
(=  m  -  q  +  1,  where  m  =*  3  and  q  ■  2)  of  the  positive  u^  variables  to  define  the 
convex  set.  Only  u^  and  Ug  are  positive  at  the  origin,  and  thus  the  convex  region 
is  uniquely  given  as  the  intersection  of  the  half  spaces  corresponding  to  (l)  and  (2). 
Extending  the  edges  from  the  origin  to  the  boundary  of  the  convex  set  gives  the  GLP  cut 
Indicated  by  the  dotted  line. 

The  new  cut  for  tills  problem  may  be  determined  by  reference  to  Remark  2,  since 
u^  and  Ug  must  be  the  starting  variables  for  the  sets  P^  and  Pg,  completing  the 
assignment  of  all  basic  u^  leaving  u^  which  is  nonbasic  at  the  LP  solution.  Thus, 
the  edge  corresponding  to  this  variable  (which  is  in  this  case  the  xg  axis)  bypasses 
the  first  hyperplane  it  intersects  (ug  *  0)  and  the  resulting  cut  is  precisely  the 
hyperplane  (l)  Itself  (i.e.,  its  associated  half  space). 

The  second  of  the  two  example  problems  is  the  same  as  the  first  with  the  additional 
constraint 

*2  <  2  • 

The  variable  u^  is  redefined  to  be  the  slack  for  this  new  constraint  (i.e.,  u^  *  2-Xg). 
The  geometric  analog  of  this  problem  appears  in  Figure  2.  The  feasible  LP  region  is 
bounded  by  (l),  (2),  (3)  (which  correspond  to  u^  ■  0,  Ug  ■  0,  u^  *  0  --  for  the  new  u^) 
and  by  the  two  coordinate  axes.  Once  again  the  origin  is  the  optimal  LP  solution,  and 
the  GLP  cut  of  [2]  specifies  the  selection  of  2  positive  u^  variables  to  determine  the 
convex  set.  This  time,  however,  there  are  three  possible  choices  since  all  three  of 
the  u^  are  positive  at  the  origin,  and  any  pair  is  a  legitimate  selection.  The  two 
best  choices  are  the  pairs  u^,  Ug  and  u^,  u^  giving  rise  to  the  cuts  indicated  by 
the  dotted  lines.  (The  third  choice,  ug#  ^  gives  the  cut  that  goes  through  the 
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intersections  of  (2)  and  (3)  with  tht  coordinate  axes,  and  hence  is  clearly  dominated 
by  the  other  cuts.) 

The  new  cut  for  this  problem  may  be  determined  by  Remark  3,  since  all  of  the  u^ 
are  basic.  The  three  choices  for  variables  to  define  the  GLF  cuts  of  [2]  provide  the 
choices  for  the  "starting"  variables  of  the  sets  and  Pg.  To  derive  the  new  cut 
we  take  P1  ■  (1)  and  P£  ■  {2},  corresponding  to  the  selection  of  the  starting  variables 
u^  and  Ug  (which  thus  become  the  initial  z^  and  Zg)-  The  basic  variable  u^  will  be 
added  to  one  of  these  sets  in  the  manner  indicated  by  Remark  Extending  the  edges 
from  the  LP  vertex,  we  note  that  the  blocking  hyperplane  (2)  (Zg  =  0)  is  intersected 
before  (3)  along  the  x^  edge,  and  thus  we  assign  u^  to  the  same  set  as  Ug,  producing 
the  completed  sets  P1  *  (1)  and  Pg  «  (2,3).  Linear  combinations  of  ug  and  u^  (to 
create  the  new  Zg)  give  the  hyperplanes  passing  through  the  intersection  of  (2)  and 
(3) •  In  this  fashion  the  first  Eg  hyperplane  (ug  ■  0)  may  be  "rotated  aside"  by  the 
addition  of  u^  to  Pg  to  produce  a  hyperplane  that  is  not  intersected  by  either  of  the 
edges  until  after  their  Intersection  with  (1).  Correspondingly,  it  is  easy  to  specify 
a  linear  combination  that  removes  Zg  from  the  category  of  a  blocking  variable  for  the 
edge  without  making  it  into  a  blocking  variable  for  the  Xg  edge.  (The  reader  may 
verify  this  by  reference  to  the  LF  tableau  for  this  example  problem.)  As  a  consequence 
both  edges  may  be  extended  to  their  intersections  with  (l),  and  the  new  cut  is  given 
by  (1)  itself.  As  in  the  previous  example,  this  is  the  strongest  cut  possible. 

It  is  interesting  to  note  that  if  the  pair  u^,  u^  were  selected  to  be  the  start¬ 
ing  variables,  then  the  same  cut  would  still  be  obtained.  (That  is,  the  strategy  of 
Remark  b  would  still  generate  the  same  sets  P^  and  Pg  and  the  same  hyperplanes 
z^  «  0  and  Zg«0  as  before  —  with  the  components  of  Pg  selected  in  reverse  order.) 
Moreover,  if  the  "worst"  pair  ug,  were  selected,  then  the  new  cuts  that  result 
by  associating  u^  with  either  of  these  variables  are  still  stronger  than  the  best  of 
the  original  GLP  cuts. 

The  increased  strength  of  the  new  cuts,  as  clearly  illustrated  by  the  preceding 
examples,  suggests  the  desirability  of  using  them  in  place  of  those  of  [2]  in  applica¬ 
tions.  _3_ 
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